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1 INTRODUCTION 

We present a simple yet rigorous theory of integration that is based on two axioms 
rather than on a construction involving Riemann sums. With several examples 
we demonstrate how to set up integrals in applications of calculus without using 
Riemann sums. In our axiomatic approach even the proof of the existence of the 
definite integral (which does use Riemann sums) becomes slightly more elegant 
than the conventional one. We also discuss an interesting connection between our 
approach and the history of calculus. The article is written for readers who teach 
calculus and its applications. It might be accessible to students under a teacher's 
supervision and suitable for senior projects on calculus, real analysis, or history of 
mathematics. 

Here is a summary of our approach. Let p : [a, 6] — > R be a continuous function 
and let / : [a, b] x [a, 6] — > E be the corresponding integral function, defined by 

I[x,y)= f p{t)dt. 

Recall that I{x,y) has the following two properties: 

(A) Additivity: I{x, y) + I{y, z) — I{x, z) for all x,y, z Cz [a, b]. 

(B) Asymptotic Property: I{x, x + h) = p{x)h + o{h) as /i — > for all x G [a, 6], 
in the sense that 

I{x,x + h)^ p{x)h ^ ^ 

h^O h 

In this article we show that properties (A) and (B) are characteristic of the 
definite integral. More precisely, we show that for a given continuous p : [a, b] — > R, 
there is no more than one function / : [a, b] x [a, 6] ^ R with properties (A) and (B). 

This will justify the simple definition p{x)dx =^ /(a, 6), where I{x, y) is a function 
satisfying (A) and (B). In this manner, we are able to rigorously develop the theory 
of definite integrals and their applications without partitioning the interval [a, b] and 
without using Riemann sums. Notice that, at this stage, the existence of the integral 
is still unsettled. Next, we prove that if R{x) is an antiderivative of then 
I{a,b) — R(b) — R{a). Conventionally, this formula is used for explicit evaluation. 
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but in our approach it also guarantees the existence of the integral for all p{x) 
with an explicit antiderivative R{x). Our approach utilizes the Riemann partition 
procedure for only one purpose: to prove the existence of the definite integral for 
functions without explicit anti-derivatives such as f{x) = e~^ . Also, our proof 
(Theorem I3.2p seems to be slightly more elegant than the conventional one. Next 
we use axioms similar to (A) and (B) to define the concepts of area under the 
curve, arclength, volumes of revolution, etc. More precisely, we define all of these 
geometrical quantities as being additive and asymptotically equal to their Euclidean 
counterparts (such as the area of a rectangle, the Euclidean distance between two 
points, the volume of a cylindrical shell, etc.). Our definitions are mathematically 
correct and well motivated. Derivations of the corresponding integral formulas (for 
the area under a curve, arclength, volume, etc.) appear in our approach as simple 
rigorous theorems with proofs done in the spirit of asymptotic analysis; we involve 
neither partitions of the interval [a, b] nor Riemann sums. 

The elementary theory of integration presented in this article (and summarized 
above) opens the door for a simple yet completely rigorous method of teaching 
integration and its applications in a calculus course or a beginning real analysis 
course. Also, we strongly recommend our method of setting up integrals (without 
Riemann sums) for teaching physics and engineering courses based on calculus. 

The method presented in this article has a long and interesting history. The 
reader might be surprised to learn that practically all textbooks on calculus and its 
applications that were written in the period between Leibniz and Riemann motivate, 
define, and set up integrals by methods very similar to the method presented in this 
article, although disguised in the language of infinitesimals. The method (among 
other treasures) was lost in the history of calculus after infinitesimals were abolished 
at the end of 19th century. This explains the choice of the title "A Lost Theorem...." 
We shall briefly discuss the connection of our approach with infinitesimal calculus 
in Section [5] In the modern literature we identify three sources using methods for 
integration similar to our approach. In H. J. Keisler's calculus textbook ^ (look 
for Infinite Sum Theorem on p. 303), the reader will find a method similar to ours 
in the framework of nonstandard analysis (in a very accessible form). We also refer 
to S. Lang (0, pp. 292-296), L. GiUman and R. McDowell ([1], pp. 176-179) and L. 
Gillman [3], where a property similar to (B) is used for a definition of the definite 
integral. 

2 TOPICS IN ASYMPTOTIC ANALYSIS 

For our axiomatic approach to integration we assume knowledge of limits, continu- 
ity, and derivatives at the level of a typical beginning calculus course. In contrast, 
we do not assume any knowledge of integration. Finally, we need several elementary 
concepts and notation borrowed from asymptotic analysis, which we present in this 
section. Most of the results are elementary, and we leave the proofs to the reader. 

Definition 2.1 We denote by o(a;") the set of all functions / : — > R such that 

Df C K, e Df, and linix^o f (x) / x"- = 0, where Df stands for the closure of 
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Df in M. We summarize this as o(x") = {f{x) : f{x)/x'^ — > as x — > 0}. It is 

customary to write f{x) = o(.t") instead of the more precise / € o(.x") in the case 
when / is an unspecified function in o{x") . If n = 0, the above definition reduces 
to fix) = o(l) if fix) ^ as X ^ 0. 

Example 2.1 x^ = o(a;) since /x ^ as x ^ 0. In contrast, since ^ o(x) since 
sina;/x — > 1 as a; — > 0. Also, we have sina; = o(l) since sina; — > as a; — > 0. 

Theorem 2.1 (Increment Theorem) Let f : [a,b] ^ R be a function and x € 
(a, b). Then f is differentiable at x if and only if fix + h) — fix) = f'ix)h + o(/i). 

Lemma 2.1 Let fix) = o(a;'") and gix) = o(x"). Then (a) fix) ± gix) = oix'^) 
where k = min(m, n), and (b) fix)gix) — o(x'""'""). Also, (c) if fix) = o(l) and 
gix) = o(l) and f is continuous at 0, then figix)) = oil). 

Remark 2.1 (Asymptotic Algebra) It is customary in asymptotic analysis to 
write simply o(x™) ± o(x") = o(x'=), o(a;'") ± o(a;") = o(x'"+") and o(o(l)) = o(l) 
instead of (a), (b), and (c) in the above lemma, respectively, when no confusion 
could arise. We shall use this notation. 

Lemma 2.2 (Absolute Value) Let AgR. Then \A + o(a;")| = \A\ + o(a;"). 

3 DEFINITE INTEGRALS WITHOUT RIEMANN 
SUMS 

Lemma 3.1 (Uniqueness) Let a,b £ R, a < b, and let p : [a,b] ^ R fee o 
continuous function. Let I : [a, b] x [a, b] ^ R be a function in two variables 
satisfying the axioms (A) and (B) at the beginning of the introduction. Then J(a, x) 
satisfies the initial value problem ^/(a,a;) = pix), /(a, a) = on the interval 
(a, 6). Consequently, there is no more than one function I : [a,b] x [a,b] — > R 
satisfying the properties (A) and (B). 

Proof: Suppose that Iix,y) satisfies (A) and (B). We have 

d Iia,x + h) - Iia,x) (A),. Iix,x + h) 
—I(a,x) = lim — i- = lim — — ; = 

dx h^o h h^o h 

(1) 

Iix,X + h)-pix)h + pix)h ^ Iix,X + h)-pix)h ^ (B) p ^ ^ 

h^O h h^Q h 

Also, (A) implies /(a, a) +/(a, 6) = /(a, 6), so /(a, a) = 0. Suppose that 
Jix,y) is another function satisfying the axioms (A) and (B), and let A(x, y) = 

lix.y) — Jix,y). Wc have ^A(a, .t) = and A(a, a) = 0, and hence A(a,a;) = 
for all X in (a, b). Consequently, for a given pix) there is no more than one /(a, x). 



3 



Next, (A) implies I{x,y) — I(a,y) — I{a,x). Thus I{a,x) uniquely determines 
/ : [a, b] X [a, b] ^ R. A 

Notice that the uniqueness result presented above does not involve partitioning 
of the interval [a, b] and Riemann sums. Rather, it is based on the more elementary 
result from calculus that "every function with zero derivative on an interval is a 
constant." 

The above lemma justifies the following definition: 

Definition 3.1 (Axiomatic Definition) Let p : [a, 6] ^ R be a continuous func- 
tion and let / : [a, b] x [a, 6] ^ M be a function in two variables satisfying the axioms 
(A) and (B) at the beginning of the introduction. Then the value I{a,b) is called 
the integral of p{x) over [a, b]. We shall use the usual notation jj^ p{x)dx '= I{a, b). 

Note that axioms (A) and (B) in the beginning of the introduction are easily 
motivated and visualized (see Figure 1). We observe as well that if a < a < /? < a, 
then the restriction / \ [a, f3] x [a,P] also satisfies the axioms (A) and (B), and 
thus p{x)dx = I{a, P). 

Theorem 3.1 (Fundamental Theorem) As before, let p{x) be a continuous func- 
tion on [a,b]. 

(i) // the integral I{a,x) of p{t) over [a,x] exists for every x in [a,b], then 
■£l{a,x) = p{x) on {a,b). 

(ii) If R{x) is an anti- derivative of p on [a, 5], then I{a,b) = R{b) — R{a). 

Proof: Part (i) follows directly from Lemma [3TT1 For (ii), we have -^I{a,x) ~ p{x) 
by (i) and^_R(x) = p{x) by assumption. It follows that -^[I{a,x) — R{x)\ = on 
{a,b), implying R{x) ~ I{a,x) + C for some constant C. Hence, R{b) — R{a) = 
[I{a, b) + C ~ {I{ci, a) + C)] ~ I{a, b), as required, since I{a, a) by Lemma [5TT1 
▲ 

Notice that our simple theory assumes that the function I(x,y), and thereby 
the integral /(a, 5), exists. The following is our first existence result (for the general 
existence result see Theorem 13. 2p . 

Corollary 3.1 (Weak Existence) If p{x) has an antiderivative on [a,b], then 
the integral of p{x) over [a, b] exists. 

Proof: Let R{x) be an antiderivative of p{x). Then the function I : [a, b] x [a, &] ^ R 
defined by I{x,y) = R{y) — R{x) satisfies (A) and (B) (at the beginning of the 
introduction) and it is the only function which satisfies (A) and (B) by Lemma l3Tl 
The number /(a, b) is the integral we are looking for. ▲ 

Remark 3.1 (Basic Properties of the Integral) The basic properties of the 
integral follow immediately from part (ii) of Theorem 13.11 under the assumption 
that the integrals exist. For example. Theorem 13.11 implies the linear property 
J[cif{x) + C2g{x)]dx — ci J f{x)dx + C2 J g{x)dx provided that at least two of the 
three integrals exist. 
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Our simple yet rigorous theory presented so far is powerful enough to support 
most of the topics related to the integral and its applications in a typical beginning 
calculus course. To deal with integrals such as dx and sin {x'^)dx, we need 

a general existence result. This (and only this) is the place in our approach where 
we use partitions of the interval [a, b] and Riemann sums. 

Let [a, b] be, as before, a closed interval in K and x,y £ R, a < x < y < b. Recall 
that a partition of [x,y] is a finite ordered set of the form P — {xg, xi, . . . , x„}, 
where n £ N and x — xq < xi < ■■■ < Xn — y. We denote by V[x,y] the 
set of all partitions of Let p : [a, 6] ^ R be a continuous function and 

P = {xq, Xi, . . . , Xn} be a partition of [a;, y]. Recall that 

L{P) = V ( niin p{t) I {xk - Xfc-i), 
C/(P) = V ( max p{t) ] {xk - Xk^i), 

\Xk-l<t<Xk J 

are called the lower and upper Darboux sums of p{t) determined by P, respectively. 
Let X, y, z e M, a < a; < y < z < fo. Notice that if P G P[a;,2/] and Q G P[y,z], 
then PUQ € P[x, z] and we have L{P) + L{Q) = L{P U Q) and U{P) + U{Q) = 
U{PUQ). The next result can be found in any contemporary textbook on Riemann 
integration. 

Lemma 3.2 Let p : [a,b] R be a continuous function. Then: 

(i) For every two partitions P and Q of [a, b] we have 

{mm p){b- a) < L{P) < U{Q) < {max p){b~ a). 

[a,b] [a, 6] 

(ii) For every e G R_|_ there exists a partition P of [a, b] such that U{P) — L{P) < 

e. 

Theorem 3.2 (General Existence Result) Let p : [a,b] ^ M. be a continuous 
function. Then p has both an integral and an antiderivative. 

Proof: Suppose, first, that x,y £ M., a < x < y < b. We observe that the set 
{L{P) I P G V[x, y]} is bounded from above by the number (max2;<t<j, p{t)){y — x) . 
Thus J(x, y) — sup{L(P) \ P € V[x,y\] exists in R by the completeness of R. We 
intend to show that /(x, y) satisfies the axioms (A) and (B) at the beginning of 
the introduction. We start with (B): we have (min2;<t<j, p(t))(?; — x) < I{x,y), by 
the definition of /(x, y), since P = {x, y} is a partition of the interval [x, y]. We let 
y~x ^ h and the result is {mmx<t<x+h p{t))h < I{x,x + h) < {mayix<t<x+h p{t))h. 
The latter implies axiom (B) since p is continuous at x. To prove (A), we observe 
that L{P) < I{x,y) < U{P) for every partition P G 7'[x,y]. Indeed, the first 
inequality follows directly from the definition of I{x, y) and the second inequality 
follows from the definition of I{x,y) and part (i) of Lemma [3.21 Next, suppose 
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that x,y,zGM.,a<x<y<z<b, and let Q G ^[2/, z]. As before we have 
HQ) < I{y, z) < U{Q), and also L{P U Q) < I{x, z) < U{P U Q) since F U Q is a 
partition of [x, z\. After summing up we obtain: 

L[P) + L{Q) - U{P U g) < I[x, y) + /(y, z) ~ I{x, z) < U{P) + U[Q) - L{P U Q). 

Now we can choose partitions P £ Vlx, y] and Q £ Vly, z] such that U{P) — L{P) < 
e/2 and U{Q)-L{Q) < e/2 by part (ii) of LemmaO Also, U{PUQ)-L{PUQ) < e 
since L{P) + L{Q) = L{PUQ) and U{P) + U{Q) = U{PUQ). Thus the above chain 
of inequalities reduces to — e < I{x,y) + I{y,z) — I{x,z) < e, implying I{x,y) + 
I{y,z) — I{x,z), as required. Finally we can eliminate the restriction on x,y and 
z by extending the function I{x, y) to a function in the form / : [a, h] x [a, 6] — > R 
by letting I{x,y) = —I{y,x) and I{x,x) = for all a; and y in [a, 6]. Notice that 
the function I(x,y) just defined also satisfies (A) and (B); thus it is uniquely 
determined, by Lemma |3. II The number /(a, 6) is the integral we are looking for, 
by Definition 13.11 and R{x) — I{a,x) is an antiderivative of p{x), by part (i) of 
Theorem 13. II ▲ 

4 SETTING UP INTEGRALS WITHOUT RIE- 
MANN SUMS 

Imagine that you are a young instructor preparing to cover arclength in a typical 
calculus course. We can safely assume that before going in front of your students 
you would like to clarify the structure of the topic for yourself: what is the definition 
of "arclength," how do I motivate it, is there a theorem to present (with or without 
formal proof), and finally which examples should I use? One option is to define the 
concept of arclength directly by the integral formula L(a, b) = -^/l + f'^{x)dx. If 
this is your choice, your next task will be to motivate this definition. You might use 
the integral formula to calculate some familiar results from high school mathemat- 
ics: the distance formula for a line segment or the formula for the circumference of 
a circle. This approach, although completely legitimate, is rarely used by calculus 
textbooks. The integral formula still looks terribly unmotivated even after deriving 
the distance formula for a line segment. Also, it is far from clear that this is the 
only formula producing the distance formula. For that reason most calculus books 
use Riemann sums to convince students that the integral formula is "reasonable." 
What follows is messy mathematics: the partition of the interval [a, b] gives the 
impression that the step of the partition Ax is "very small" which leads to the 
conclusion that AL is approximately equal to Ax'^ + Aj/^. The latter produces 
our integral formula after factoring out Ax and taking the limit as Ax — > 0. You 
might ponder for hours questions such as: a) What, after all, is the definition of 
"arclength"? Is the integral formula exact or approximate? After all, its derivation 
involves the approximate formula AL « y^Ax^ + Ay^. b) What is meant by "Ax 
is very small"? If Ax = 0, then Ay — and the root Ax"^ + Ay"^ is also equal to 
zero. If Ax ^ 0, then AL ^ ^ Ax^ + Aj/^ (unless the curve is a straight line), c) Is 
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the "derivation" of the formula for arclength a sort of casually presented proof of a 
casually stated theorem? And if yes, what is the rigorous version of this theorem? 
Worst of all is that the students are usually preoccupied with the technicalities of 
the partition procedure and the sigma notation in the Riemann sums and hardly 
pay attention to the fact that a new important concept has just been introduced. 

In this section we take another approach. The concept of arclength is defined 
as an additive quantity which is asymptotically equal to the Euclidean distance 
between two points (Definition 14. 2p . The definition is mathematically correct and 
well motivated. It is based on the concept of limit, not integral. Next we derive the 
arclength integral formula as a simple rigorous theorem in the spirit of asymptotic 
analysis (Theorem I4.2p . Similarly we define the rest of the additive quantities from 
geometry and physics such as area under a curve, volume of a solid of revolution, 
etc. We involve neither a partition of the interval [a, h] nor Riemann sums. 

We start with area under a curve. Our assumption is that the reader knows 
what the area of a rectangle is but not necessarily what the area under a curve 
is. In particular, we do not assume that the reader necessarily knows the integral 
formula A{a, b) — f{x)dx for the area under a curve; rather our goal is to derive 
this formula starting from the more elementary concept of the area of a rectangle. 

Definition 4.1 (Area Under a Curve) Let y — f{x) be a continuous function 
on [a, b] such that f{x) > for all x S [a, h]. Let A : [a, b] x [a, 5] ^ R be a function 
in two variables satisfying the following two properties: 

(a) A{x, y) + A{y, z) ~ A(x, z) for all x,y, z £ [a, b] . 

(b) A{x, x + h) — ±R{x, x + h) + o(h) as h ^ Q± for all x E [a, 6], where R{x, x + h) 

denotes the area of the rectangle with vertices (x, 0), {x + h, 0), {x + h, f{x)), 
and [x, f{x)). 

The number A{a, b) is called (by definition) the area under the curve y ~ f{x) 
and above the interval [a, b]. 

The above definition can be easily motivated (see Figure 2). 

In the next theorem we derive the familiar integral formula for A{a, b) without 
partitions or Riemann sums. While deriving this formula, we shall simultaneously 
prove two things: (a) the correctness of the above definition, and (b) the existence 
of the area under the curve. As in the conventional approach, the integral formula 
offers a practical method for explicit evaluation. 

Theorem 4.1 A(a,b) = ^^f{x)dx. Consequently, the area A(a,b) is uniquely 
determined by the properties (a) and (b) in Definition \4^.1\ 

Proof: We have to find the asymptotic expansion of A{x, x + h) in powers of h as 
h ^ and extract the coefficient in front of h. Since R{x,x + h) = f{x)\h\, we 
have A(x,x + h) = ±R{x,x + h) + o{h) = ±f{x)\h\ +o{h) = f{x)h + o{h), and the 
above formula follows directly from Definition 13.11 for p{x) — f{x). Conversely, it is 
easy to verify that the function A : [a, b] x [a, 6] ^ R defined by A{x, y) — f{t)dt 



7 



satisfies the properties (a) and (b) in the above definition, thus the number A{a, h) 
is the area under the curve. ▲ 

Next, we define the concept of arclength without partitions or Riemann sums. 
Our assumption is that the reader knows what the Euclidean distance between two 
points is but not necessarily what the arclength along a curve is. In particular, we do 
not assume any knowledge about the integral formula L(a, b) — J^^ -y/l + f'^{x)dx; 
our goal is to derive this formula starting from the more elementary concept of 
Euclidean distance between two points. 

Definition 4.2 (Arclengtii) Let / e C^[a,b] and let L : [a,b] x [a, 6] R be a 
function in two variables satisfying the following two properties: 

(a) L{x, y) + z) = L{x, z) for all x,y, z Cz [a, b]. 

(b) L{x, x + h) = ±D{x, x + h)+o{h) as /i ^ 0± for all x € [a, 6], where D{x, x + h) 

is the Euclidean distance between the points (a:, f{x)) and {x + h, f{x + h)). 

The number L{a, b) is called (by definition) the arclength of the curve y = f{x) 
between the points (a, /(a)) and {b,f{b)). 

The above definition can be easily motivated (see Figure 3). 

In the next theorem we rigorously derive the formula for arclength without 
involving partitions of the interval or Riemann sums. Among other things we prove 
correctness of the above definition and the existence of the arclength L(a, 6). 



Theorem 4.2 L{a,b) — 1^1 + f'^{x)dx. Consequently, the arc length L {a, b) is 
uniquely determined by the properties (a) and (b) in Definition \4-S\ 

Proof: We have to find the asymptotic expansion of L(x,x + h) in powers of h 
as ft, and extract the coefficient p{x) in front of h (see Definition 13. ip . Let 
Ay = f{x + h) — f{x) and recall that Ay = f {x)h + o{h) (Increment Theorem l2.ip . 
We have: 



L{x, x + h)= ±D{x, x + h) + o{h) = ±\/ h'^ + Ay"^ + o{h) 



±W,;i + |'^'W" + °""V + o(,.) 



= h^i + r^{x)+o{i) + o{h) = h [vi + /"(•^) + + /"(^) + o(i) - v^tthx) 

0(1) 



+ o{h) = h 



^i + nx) + 



v/l + /'2(a;)+o(l) + v/l + /'2(x) 



+ o{h) = 



= h 



v/1 + /'2(.t)+o(1) + o{h) = hy/T+fHx) + o{h) + o{h) 



h^l + THxj + o{h). 
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For the coefficient in front of h we fiave p{x) = ^^1 + /'^(a;), which implies our 
integral formula by Definition 13. II Conversely, it is easy to verify that the function 
L : [a, b] x [a, 6] — > R defined by L{x, y) = /J i/l + f''^{t)dt satisfies (a) and (b) in 
Definition 14.21 . Thus the number L(a, b) is the arc length of the curve. ▲ 

Next, we set up the integral for a volume of revolution about the y-axis without 
partitions of the interval and Riemann sums. Our assumption is that the reader 
knows what a volume of a cylindrical shell is but not necessarily what a volume of 
revolution is. 

Definition 4.3 (Volume of Revolution) Let y — f{x) be continuous on [a, 6], 
fix) > 0, and < a < 6. Let V : [a.,b] x [a,b] — > M be a function in two variables 
satisfying the following properties: 

(a) V{x, y) + V{y,z) = V{x, z) for all x, y , z e [a,b]. 

(b) V{x,x + h) = ±U{x,x + h) + o{h) as h ^ 0± for all x S [a,b], where 
U{x,x + h) is the volume of the cyhndrical shell obtained by revolving the rectangle 
with vertices (x, 0), {x+h, 0), {x+h, f{x)) and [x, f{x)) about the y-axis (see Figure 
4). 

The number V{a, b) is called (by definition) the volume of revolution about 
the y-axis of the curve y — f{x). 

In the next theorem we rigorously derive the familiar formula for the volume 
V{a^ b), and we do so without partitions or Riemann sums. Among other things we 
prove the uniqueness and existence of the volume V{a^ b). 

Theorem 4.3 V{a, b) = 2'Kxf{x)dx. Consequently, the volume V{a, b) is uniquely 
determined by the properties (a) and (b) from Definition \^.3\ 

Proof: The volume of the cylindrical shell is x + h) ~ \tt{x + /i)^ — 7rx^|/(x). 
Hence, with the help of Lemma [2.21 we have V{x, x + h) = ±U (x, x + h) + o{h) = 
±\tt{x + hf - 7ra;2|/(x) + o{h) = ±\n{2xh + h'^)\f{x) + o{h) = ±2TTxf(x)\h\ + 
TTf{x)h'^ + o{h) = ±2TTxf{x)\h\ + o{h) + o{h) = 27rxf{x)h + o{h), and from Defini- 
tion [SH] we have p{x) = 2TTxf{x). Conversely, it is easy to verify that the function 
V : [a,b] X [a,b] R defined by V{x,y) = 2ntf(t)dt satisfies the properties (a) 
and (b) in Definition 14.31 Thus the number V{a,b) is the volume of revolution. A 

5 BACK TO INFINITESIMAL CALCULUS: A 
LOST THEOREM 

As we explained at the end of the introduction, the method presented in this article 
has a long and interesting history. The purpose of this section is to establish a rela- 
tion between our method of integration and infinitesimal calculus. This section may 
be helpful to those readers who are interested in reading original texts on infinitesi- 
mal calculus but who might not have background in modern nonstandard analysis. 
We should mention that the tradition of using infinitesimal arguments is still very 
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much alive and can be found in many contemporary texts on applied mathematics, 
physics, and engineering. For example, in the famous The Feynman Lectures on 
Physics [T] we located about ten cases of integral formulas for additive physical 
quantities derived in the spirit of infinitesimal calculus and without Riemann sums 
(see Volume I: pp. 13-3, 14-8, 43-2, 44-10 and 44-11, 46-6, 47-5; Volume II: pp. 3-2, 
38-6). For that reason we believe that our article and the discussion in this section 
in particular might also be helpful to pure mathematicians who are interested in 
reading texts on applied mathematics, physics, or engineering but who might feel 
uneasy with infinitesimal reasoning. 

Recall that in the period from Leibniz to Weierstrass, calculus was commonly 
known as infinitesimal calculus and was based on the hypothesis that there exist 
nonzero infinitesimals, i.e., mysterious numbers dx with the property < \dx\ < 1/n 
for all n e N. We should keep in mind that at that period not only the theory of 
infinitesimals but also the theory of real numbers was without rigorous foundation. 
So, the existence of nonzero infinitesimals should not be dismissed as nonsense. Yes, 
the field of the real numbers M does not have nonzero infinitesimals, but there were 
no real numbers in the era of Leibniz and Euler; the real numbers were an invention 
of the late 19th century and were systematically implemented in mathematics at 
the beginning of 20th century. 

We shall discuss Leibniz-Euler infinitesimal calculus using the common "differ- 
ential notation" : Let y = f{x), a < x < b he a, function. In what follows, dx stands 
for a new independent variable (real or infinitesimal depending on the context) and 
Ay — f{x + dx) — f(x) stands for the increment of y. If / is differentiable at x, 
then dy = f'{x)dx stands for the differential of y. It is clear that dx = Ax, but we 
prefer to use dx over Ax, thus keeping track of the fact that x is an independent 
and y is a dependent variable. Before proceeding further we should notice that in 
the old infinitesimal calculus, and in the modern nonstandard analysis as well, the 
notation dx rarely stands for a fixed infinitesimal number (as tt, e, ^/2, etc. stand for 
specific irrational numbers). Rather, dx is usually used for an independent variable 
ranging over a set consisting of both infinitesimal and real (standard) numbers. To 
demonstrate this point we shall present the characterization of continuity used by 
Euler but rigorously justified in the modern nonstandard analysis: Let a; be a real 
(standard) number in the domain of /. Then / is continuous at x if and only if 
f{x + dx) — f{x) is infinitesimal for every infinitesimal dx. In the manuscripts of 
Euler this statement appears in a slightly more casual form: / is continuous at x if 
and only if Ay is infinitesimal whenever dx is infinitesimal. It is clear that if both 
X and dx are real (standard) numbers, then Ay — f{x + dx) — f{x) is also a real 
(standard) number. Also, "real (standard) number" is a modern term. Leibniz and 
Euler would rather use "usual quantity" (as opposed to "infinitesimal quantity" ) 
instead. 

Let us try to mimic, for example, the arguments used by Euler for setting 
up the integral for the arclength L{a,b) along the curve y = f{x),a < x < b 
in a typical calculus course in the middle of the 18th century (compare with our 
Definition 14.21 and Theorem 14. 2p . If dx is infinitesimal, the arclength L{x,x + dx) 
between the points {x, f{x)) and (x+dx, f{x + dx)) should be equal to the Euclidean 
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distance D{x, x + dx) — \/ dx-^ + Ay^ between the same points, up to infinitesimals 
of second order relative to dx; in symbols, L{x,x + dx) « dx^ + Ay^ (see Figure 
3, where h should be replaced by dx). On the other hand, since L is an additive 
quantity, we have L(x, x + dx) = L(a, x + dx) — L{a, x). Thus L(x, x + dx) ~ dL, 
by the increment theorem (in its original infinitesimal form). Also, Ay « dy, by 
the increment theorem. The result is dL « -^/l + (dy/dx)'^ dx which implies the 

familiar integral formula L{a,b) — jj^ -^/l + (dy/dx)^ dx. We should note that the 
implication in the previous sentence has never been rigorously justified in the old 
infinitesimal calculus; among other things the goal of our article is to fill this gap. 

A contemporary mathematician, unless familiar with nonstandard analysis, will 
have difficulty recognizing our asymptotic property (B) (at the beginning of the 
introduction) using Euler's language presented above. For those who are interested 
in using nonstandard analysis we recommend H. J. Keisler 4 or T. Todorov [Q, 
where the reader will find more literature on the subject. In what follows, however, 
we shall choose another path: we shall use the language of asymptotic analysis 
(Section [2|) to relate the method of integration presented here with the method of 
infinitesimal calculus. For that purpose we suggest the following modification of 
Section [3] and Section |4] of our article: 

1) First, instead of the letter h (used in our article so far) we shall use the 
original Leibniz notation dx. In other words we let dx = h and we treat dx as a 
new real (standard) independent variable. In this notation our Asymptotic Property 
(B) (see the beginning of the introduction) becomes 

(B') I{x, X + dx) — f{x)dx + o{dx) as dx for all x £ [a, b]. 

As before, f{dx) = o{dx^) means f{dx)/dx" ^ as ^ 0. 

2) Our next goal is to give precise meaning to the relation w used by Euler in 
our previous example. Let F{dx) and G{dx) be two real functions. We say that 
F{dx) and G(dx) are equal up to infinitesimals of second order relative to dx, in 
symbols, F w G, if F{dx) - G{dx) = o{dx) (i.e., if {F{dx) - G{dx))/dx ^ as 
dx ^ 0). The property (B') becomes 

(B") I{x,x + dx) ~ f{x)dx for all X G [a, 6]. 

3) We might stop writing dx — > 0, since the symbol dx is more than suggestive. 

4) Now we can rewrite Section[31 replacing h by dx and the axiom (B) by (B"). 
The additive property (A) (at the beginning of the introduction) does not need 
modification. 

5) Finally, we have to replace all parts (b) in the definitions in Section [4] by 
their (b")-counterparts in the spirit of (B"). For example, part (b) of Definition l4.2l 
should be replaced by: 

(b") L{x, X + dx) w zLD{x, x + dx) for all x G [a, b], 

and similarly for the rest of the (b)s in Section 21 

While preserving the content of the article, the modification suggested above 
makes the arguments of the old infinitesimal calculus more transparent. 
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6) Readers who feel uncomfortable with the notation dx (and especially with 
dx 0) should restore the notation h used in Section 3 and Section 4 of this 
article. However we recommend our trick "replace h by dx" to those readers who 
are interested in reading original texts on infinitesimal calculus, but do not have 
background in the modern nonstandard analysis. 

Among other things our article suggests that the reasoning of mathematicians 
in the era of Leibniz and Euler, as well as the reasoning of contemporary applied 
mathematicians and physicists, is often more reliable and rigorous than is usually 
acknowledged by pure mathematicians. 

Acknowledgement: We are thankful to our colleague Donald Hartig who 
made several useful remarks on the manuscript. 
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